Distortion blockade in classical nano-electromechanical resonator 
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We consider a single electron transistor where the central island can oscillate. It has been shown 
that for weak coupling of the elastic and electric degrees of freedom the position of the island fluc- 
tuates with a small variation of the current through the device. In this paper we consider the strong 
coupling limit. We show that the system undergoes a static mechanical instability that is respon- 
sible for the opening of a gap in the current voltage characteristics even at the degeneracy point. 
We provide an analytical description of the transition point. We also discuss how the mechanical 
nature of the suppression of the current can be probed experimentally by a slow modulation of the 
gate voltage. 

PACS numbers: 



Nanoelectromechanics constitutes a rapidly developing 
and promising field of mesoscopic physics [H, @, Q ■ A par- 
ticularly important and investigated device is the single 
electron transistor with mobile parts 0, H, @|. Due to 
the Coulomb blockade the interplay between the elec- 
trical and mechanical degrees of freedom influences the 
current in a sizable way. When the oscillation ampli- 
tude is sufficiently large with respect to the tunnelling 
length, the modulation of the tunnelling rates leads to 
the shuttle phenomena, where the oscillations becomes 
synchronized with the electron tunnelling 0, H, 0, Hcj |. 
If the oscillation does not modify the distance between 
source and drain, it may nevertheless induce a modula- 
tion of the tunnelling rates through the variation of the 
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electric potential energy in space 
Typically this effect is due to the position dependence of 
the gate capacitance. If the Coulomb force generated by 
the variation of the number of electrons by one is F, the 
distance of the two equilibrium positions of the island 
will be X = F/k where k is the effective elastic constant 
of the island. The corresponding variation of the elastic 
energy is Ee = F 2 /k. For source-drain bias voltages V 
of the order of EeI& (with e the electron charge), the in- 
terplay between the electrical and mechanical degrees of 
freedom becomes important. Typically Ee is very small, 
but within the range of observation. For a nanotube of 
500 nm length, 1 nm radius, suspended at 100 nm from 
a gate Ee is of the order of 100 [ieV . Previous analytic 
work concentrate on the case eV » Ee [H, [Hj|. Very 
recently the current for eV > Ee at the degeneracy point 
has been considered numerically [Hj]. In this paper we 
present an analytic theory that gives a good description 
of the whole region eV > Ee and arbitrary gate voltage. 
We find that sweeping the gate voltage for eV > Ee the 
current jumps from in the blocked regions to a finite 
value in the conducting regions. We also discuss the be- 
havior of the system for eV < Ee- We argue how a static 
mechanical instability blocks the current for any value of 
the gate voltage. 





FIG. 1: Left panel: Circuit diagram for the single electron 
transistor with an oscillating central island. The variable X 
parameterizes an effective position. Right panel: conducting 
regions shaded in the V g -V near a degeneracy point between 
N and N + 1 electrons. 



Let us consider a single electron transistor with a mo- 
bile central island (see Fig. [1]). The model has been 
discussed in details in previous works (see for instance 



Refs. 



12, [HI, Q, [HI), we discuss the simplest con- 
figuration of symmetric capacitances, resistances, and 
bias at zero temperature. (Strictly speaking we assume 
Ee 3> ksT ^> Hlu.) For small oscillations gate capaci- 
tance varies linearly in the displacement of the island X . 
Charge transport can be described as usual by a master 
equation where the rates are modulated by the position 
of the island. For positive source drain voltage (V > 0), 
and the gate voltage V g tuned close the degeneracy point 
between TV and N + 1 electrons on the island, the non 
vanishing rates read: 

Tl.r = f(eV/2±/3 g (0)eV g T(2N+l)E c (0)TFX) , (1) 

where and Tr are the rates for transfer of electrons 
from the left lead to the island and from the island to 
the right lead, respectively, Rt is the tunnelling resis- 
tance of the left and right contacts, Ec(X) = e 2 /2Cs is 
the Coulomb energy of the device as a function of the po- 
sition of the island X, C and C g (X) are the junction and 
gate capacitances, respectively, C^(X) = 2C + C g {X), 
f(E) = E/(e 2 R T ) for E > and for E < 0, and finally 
j3 g (X) = C g (X)/C s (X). Choosing X = as the equilib- 
rium position for the island with N electrons, the force 
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acting on the island when an additional electron is added 
is then F = -(2N+l)E c [l+2C/C g (0)]C£ 1 dC g (X)/dX, 
where we set eV g = (2N + l)E c /(3 g , since we will con- 
sider small variation of V g around this value. The motion 
of the island is described by the Newton equation: 



mX(t) = -kX(t) - Fn(t) 



(2) 



where n(t) fluctuates stochastically between and 1 ac- 
cording to the rates given in Eq. (fT]) ; m is the effective 
mass of the island. It has been shown that the coupling 
to the electronic degrees of freedom introduces an intrin- 
sic damping coefficient 12j . We thus do not introduce an 
extrinsic dissipation into Eq. ([2]), assuming that the in- 
trinsic dissipation is dominant. It is convenient to intro- 
duce reduced variables: x{r) = X/X Q , u = X/(u> X ), 
v = eV/E E , v g = {P g eV g - (2N + l)E c )/E E , r = w„t, 
r o = E E /(e 2 R T w ), with uj 2 = k/m. For fixed x the 
stationary current is given by the usual expression 
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L-L r 



v 2 /4 - (v g - a;) 2 



r t v r t v {r L + r R ) 



(3) 



In the plane v-v g it leads to the Coulomb diamond struc- 
ture of the current, as shown in Fig. [1] with the degen- 
eracy point sitting at v g = for x = 0. 

The effect of the oscillation of the central island on the 
current can be studied analytically far from the bound- 
aries of the conducting regions [12j. In that limit it has 
been shown that the probability distribution of x and u 
becomes gaussian with a width controlled by the volt- 
age bias v, and the variation of the current due to the 
mechanical coupling is always a small part of the unper- 
turbed value. In this paper we will focus on the regions 
near the degeneracy point and near the two lines of tran- 
sition from the conducting to the non conducting region 
(see Fig. [TJ right panel). We assume that r o S> 1 so that 
the problem can be tackled by exploiting the separation 
of time scales between the slow mechanical oscillations 
and the frequent electronic hopping. 

Wc begin with a very simple description of the stochas- 
tic force by substituting into Eq. @ the average occupa- 
tion number: n(t) — > n(x), where: 



n(x) 



v/2 + v g — x 



for \v g —x\ < v/2 , (4) 



n(x) = for v g —x < —v/2, and n(x) = 1 for v g —x > v/2. 
This induces an average force on the island that depends 
on the occupation of the island itself. It is convenient to 
introduce an effective potential to describe this force (in 
units of Eq): 



Ueff(x) 



i(x)]dx 



(5) 



where x m = (v/2 + v g )/(l — v). The form of n[x) given 
above implies that d 2 U e ff/dx 2 equals (1—l/v) for \v g — 




FIG. 2: Left panel: Effective potential U e ff for v g = —1/2 
and v = 1.05, 1, and 0.95. Right panel: Shape of the conduct- 
ing regions (grey) in the plane v a -v when the displacement of 
the mobile part is taken into account. The small dashed dia- 
mond indicates the mechanically bistable region. 



x\ < v/2 and 1 otherwise. For v > 1 the potential has 
thus a single minimum at x m with — 1 < x m < and (by 
construction) U e ff(x m ) = 0. For < v < 1 instead two 
minima may be present, depending on the value of v g . 
The evolution of the potential as a function of v is shown 
in Fig. [5] for v g = —1/2. At v = 1 the potential becomes 
flat and for v < 1 it develops two side minima. 

The form of the effective potential has important con- 
sequences on the current flowing through the device. 
Within this simple approximation the average position 
x is determined by the minimum of the potential. The 
current is then given by Eq. with x = x m . Note that 
Xm < 0, this implies that the displacement increases the 
gate voltage seen by electrons [cf. Eq. (JXJ) ] . Thus if for 
v > 1 we sweep v g from negative to positive values, n = 
for v g < —v/2, n then starts to increase linearly from 
to 1. Due to the displacement of the average position n 
equals 1/2 (degeneracy point) at v g = —1/2 instead of 
v g =0. For the same reason it becomes \dXv g =v/2 — 1 
instead of v g = v/2. The shape of the conducting re- 
gion changes thus slightly for v ^S> 1, as shown in Fig. [5] 
But for v < 1 the consequence is more dramatic. The 
evolution in v g is no more smooth, due to the presence 
of two minima, the system is hysteretic and the stable 
position and occupation jump from x rn = — 1 and n = 
to x m — and n = 1, or viceversa. The current is thus 
always blocked for v < 1. As shown in Fig. [2] right panel, 
a gap appears in the v g -v plane at the degeneracy point. 
This effect can be regarded as the classical counterpart 
of the Frank-Condon effect in the quantum limit that has 
been shown to lead to a current suppression in molecular 
devices [l7, 3 1|1 2(|. The fact that this phenomenon 
exists also in the classical limit, makes it probably more 
easy to be observed in finite temperature device. More- 
over, the simplicity of the classical description affords 
an analytical solution and the the possibility of treating 
more complex situation, as the non stationary evolution 
discussed in the following. 

Let us now discuss the behavior of the system in the 



FIG. 3: Lines of equal energy and thus of equal probability 
Q in the x-u plane for the values of v and v g indicated. 
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conducting region (v > 1). The dependence of the cur- 
rent and the probability distribution for x along the de- 
generacy line (v g = —1/2, and v > 1) have been consid- 
ered very recently in Ref. [lbI ] numerically. We present 
here an analytical theory that gives a good description 
of the current in the whole conducting region, including 
near the boundaries and at the apex. For r o 3> 1, the 
electrons have the time to hop in and out of the grain 
many times before it can move of a sizable distance. One 
can then write a Fokker-Plank equation for the probabil- 
ity Q(x,u) [H: 

f -^^(^-^(XHQ-A.Q+^^Q (6) 

where F e (x) = —dU e ff/dx, and rj.i(x) = (dn/dv g )/vT 
is the intrinsic damping [l2l , [l4| . The non vanishing sec- 
ond moment is S(x) = J dr (<5ri(r)i5n(r)) = 2n(x)(l — 
n(x))/T Q v. It is convenient to simplify further the equa- 
tion by exploiting the fact that the fluctuating and dissi- 
pative part is small for r o 3> 1. The distribution function 
is thus a function of the effective energy: 

E e (x,u) = U eff (x)+u 2 /2. (7) 

From Eq. ([6]) we can derive an equation for V(E) = 
J dxduQ(x,u)5(E — E e (x,u))\ 



dV 
~dt 



d_ 
~dE 



-a(E)V(E) + A (P(E)V(E)) 



(8) 



with a(E) = {S(x)/2-t 1i (x)il 2 ) e and 0(E) = 
(5'(x)m 2 ) £; /2, where the averages are taken on the tra- 
jectories in the x-u plane at fixed energy E. These tra- 
jectories are shown for some values of v and v g in Fig. [3] 
Near the apex and the borders the shape is quite dif- 
ferent from an ellipse, which is the harmonic oscillator 
trajectory. 

The stationary solution of Eq. © is 



FIG. 4: Current as a function of v at v g — —1/2 from the 
analytical calculation (full line), the weak coupling theory 
(dashed line) and the Monte Carlo simulation (points). The 
values of r o are 10 (square), 1 (triangle), 0.1 (cross) and 0.01 
(full square). In the inset: current as a function of v g for 
v = 1.5. The long dashed line is Eq. (|3} for x = x m . 



where J\f is a normalization constant. Once V(E) is 
known, all physical quantities can be calculated by av- 
eraging first over the trajectories, and then over the en- 
ergy. Let us test the theory at the strongest value of 
the interaction, at the apex of the new Coulomb dia- 
mond: v = 1 and v g = —1/2. The orbits are quite simple 
for this case (see Fig. The velocity is constant for 
\x + 1/2| < 1/2 and the trajectory is exactly half an el- 
lipse for I x + 1/2 1 > 1/2. The averages of rji and S on 
this trajectories lead to the following expressions for the 
probability 



V{E) = Af(l + nV2E)e~ 6E /VE. 



(10) 



The current at the apex can be obtained by integrating 
expression (J3J with the distribution (fTD"]) . We find that 
R T I(v = l,v g = -l/2)/V = 1/(6 + 2V37r) w 0.0837. 
With the same technique one can calculate the rms of x: 
(40 + 3 V6 7r 3 / 2 ) 1/2 /(96 + 8 V6tt 3/2 ) 1/2 « 



x) 2 ) 



V{E) =AAexp 



dE'a(E')/(3(E') /(3(E) (9) 



0.628. 

The method can be applied to describe the transi- 
tion from the blocked to the conducting state at finite 
v. In order to keep the calculation simple we consider 
v 1. For v g = —v/2 the minimum of the potential 
is at x = 0. The potential is harmonic for x < with 
oscillation period 27r and harmonic for x > with period 
27r(l — 1/v)^ 1 / 2 . The oscillations will thus be more elon- 
gated inside the conducting regions (see Fig. [3]). The 
calculation at leading order in l/v gives the distribu- 
tion: V(E) =AA[l-37rVS/(8wV2(l-l/w) 1/2 )] 4(1 - 1/t,) ". 
With this distribution we find that the current jumps 
at v g = —v/2 (and symmetrically at v g = —v/2 — 1) 
with R T I(v,v g = -v/2+)/V = (£ B /eV)8/(37r 2 )[l - 
l/2v — 1/Av 2 + o(l/v 3 )]. For arbitrary values of v g and 
v we found the analytical expressions for a and (3 and 
they have been used to obtain the current by numer- 
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ical integration. To check these results we have per- 
formed a Monte Carlo simulation of the model (see for 



o 



instance Refs. [10|, |15[). Fig. 2] shows the comparison 



for v g = —1/2 and v varied between and 8, and for 
fixed v = 1.5 and v g varied between the blocked region 
to the degeneracy point (inset). We considered values 
r o ranging from 10~ 2 to 10. Surprisingly, the agreement 
between the analytical and the numerical calculations is 
very good also for r o <C 1. Apparently the current de- 
pends very weakly on r o (apart a trivial linear scaling) . 

Concerning the v g dependence of the current (inset of 
Fig. [3} compare the exact result to the weak coupling re- 
sult of Ref. 12j extrapolated to strong coupling (dashed 
line) or to the expression (|3|) taken at x = x m (Ion 
dashed line). For v g > —v/2 the current flows through 
the device inducing a fluctuation of the position. The 
observed current is then the average of Eq. © over the 
values of x visited by the island. This average increases 
the current near the threshold, producing the disconti- 
nuity, and reduces it near the degeneracy point. Also 
for the rms of the position the analytical estimate given 
above agrees with the numerical simulation: for r o vary- 
ing in the same range of Fig. 5] we find that y/ (x — x) 
ranges between 0.59 and 0.64, that compares well with 
the analytical result 0.628. The good agreement of the 
analytical and Monte Carlo results indicate that the an- 
alytical picture is accurate. It thus provide a simple and 
faithful description of device dynamics. 

The presence of the mechanical bistability can be ex- 
perimentally probed for v < 1 by modulating the gate 
voltage around -1/2: v g = — 1/2 + v° sm(u>t). According 
to the picture given before if v° < (1 — v)/2, no cur- 
rent should flow, since the island remains in the same 
stable or meta-stable position all the time: v g varies in- 
side the dashed diamond of Fig. [2j For v° > (1 — v)/2 
the island is instead always released from the metastable 
state, and it is free to oscillate between the two positions 
x = — 1 and x = 0. This implies that some current can 
flow through the device. The resulting DC current as a 
function of v° is thus discontinuous at v° = (1 — v)/2. 
In order to verify this idea quantitatively we resort to a 
Monte Carlo simulation for this non stationary situation. 
Results are shown in Fig. [5] for v = 1/2 and lo/uj = 0.1. 
For v° 3> (1 — v)/2 = 1/4 we find that the current is 
one electron per cycle, regardless of the value of r o . The 
value of the discontinuity depends instead strongly on r o . 
For r o 3> 1 many electrons can flow during the swing of 
the oscillator between the two stable solutions. In the 
opposite limit of T a <C 1, the discontinuity is strongly 
reduced. Nevertheless also in this case the nanomechan- 
ical nature of the current can be probed by studying the 
DC current as a function of the external frequency lo/lj . 
The inset shows the frequency dependence of the current 
at the threshold (v° = 1.05/4) for large and small r o . In 
both cases dips are present when lo/lj = 0.5, 1 or 2, in- 
dicating the resonant response of the mechanical degree 




g FIG. 5: DC current response to an oscillating gate voltage 



-l/2+i;° sin(ujt) as a function of v°. 



1/2,uj/lj = .1, 



(same symbols of Fig. [5} The arrow indicates the threshold 
to the conducting region. Inset: Frequency dependence of the 
DC current for v° — 1.05 the threshold value. Structures at 
lj/uj = 0.5, 1 and 2 are clearly visible. 



of freedom. 

In conclusion we have shown that the coupling to a 
mechanical degree of freedom can lead to a current sup- 
pression at the degeneracy point. A bistability appears 
for eV < Ee and we have shown how this can be de- 
tected by using an AC gate voltage. We also presented 
an analytic theory that allows to calculate with quanti- 
tative accuracy the current in the whole conducting re- 
gion (eV > Ee)- A gap at the degeneracy point has 
been already observed exp erimentally in different nano- 
mechanical devices 0, [H|. The miniaturization of the 
devices will allow to increase further the value of Ee- 
Our predictions can be useful to test if the observed phe- 
nomena are really related to nanomechanical effects. 

We acknowledge support from the French Agence Na- 
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(Grenoble). 
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